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Observation of B Meson Decays to ljK* and Improved Measurements for cup and u>fo 
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Wo present measurements of B meson decays to the final states uiK* , top, and cjfo, where K* 
indicates a spin 0, 1, or 2 strange meson. The data sample corresponds to 465x10® BB pairs 
collected with the BABAR detector at the PEP-II e+e^ collider at SLAC. B meson decays involving 
vector-scalar, vector-vector, and vector-tensor final states are analyzed; the latter two shed new 
light on the polarization of these final states. We measure the branching fractions for nine of these 
decays; five are observed for the first time. For most decays we also measure the charge asymmetry 
and, where relevant, the longitudinal polarization /l. 



PACS numbers; 13.25.Hw, 12.15.Hh, ll.30.Er 



4 



Studies of vector- vector (VV) final states in B decays 
resulted in the surprising observation that the longitudi- 
nal polarization fraction in _B — > (l)K* decays is ^-^0.5, 
not '^l The latter value is expected from simple 
helicity arguments and has been confirmed in the tree- 
dominated ^ B ^ pp decays [§| and i?+ wp+ decays 
01 . It appears that the expectation is correct for 

tree-dominated decays but is not generally true for decays 
where & — *■ s loop (penguin) amplitudes are dominant. 

There have been numerous attempts to understand the 
polarization puzzle (small /l) within the Standard Model 
(SM) [5i], and many papers have suggested non-SM ex- 
planations 0. The SM picture improved recently with 
the calculation of for most charmless VV decays [2;] 
with inclusion of non-factorizable effects and penguin an- 
nihilation amplitudes. Improved understanding of these 
effects can come from branching fraction and /l measure- 
ments in decays such as i? — > toK*, which is related to 
B (jjK* via SU(3) symmetry Among these decays, 
there is evidence for only B^ —+ loK*^ [1, . Information 
on these and related charmless B decays can be used to 
provide constraints on the CKM angles a, /3, and 7 Q. 

Further information on the polarization puzzle can 
come from measurements that include the tensor meson 
-ftr|(1430). A measurement of the vector-tensor {VT) de- 
cay B 0/^1(1430) [10'] finds a value of /l inconsistent 
with 0.5 (but consistent with 1), so a measurement of 
the related decay B cj/'ir|(1430) would be interesting. 
The only theoretical predictions for these modes are from 
generalized factorization calculations [11]; the branching 
fraction predictions for the B (jjiir|(1430) decays are 
~ (1 - 2) X 10"^, but there are no predictions for f^. 
There have been a variety of measurements for similar B 
decays that include the scalar meson _firQ(1430) [lol. IT3|. 
For the scalar- vector (SV) decays B a;ifn(1430), there 
are recent QCD factorization calculations [l3| that pre- 
dict branching fractions of about 10^^. 

We report measurements of B decays to the final states 
wK*, Lup, and w/o(980), where K* includes the spin 
0, 1, and 2 states, X^(1430), i^*(892), and K^iUSO), 
respectively. While a complete angular analysis of the 
VV and VT decays would determine helicity ampli- 
tudes fully, because of the small signal samples we mea- 
sure only /i. Given our uniform azimuthal acceptance, 
we obtain, after integration, the angular distributions 
d'^T /{d cos 9 id cos 62) ■■ 



/t sin 6*1 sin 6*2 + 4/lcos 6'iCos^2, 



(1) 



/t sin^ 6*1 sin^ 62 cos^ 6*2 + ^ cos^ 9i (3 cos^ 6I2 - 1)^ (2) 

for the VV and VT [l^l decays, respectively, where 
fr — ^ ~ Jl and di and 62 are the helicity angles in the 
V 01 T rest frame with respect to the boost axis from the 
B rest frame. For decays with significant signals, we also 
measure the direct CP-violating, time-integrated charge 



TABLE I: Selection requirements on the invariant masses and 
helicity angles of B-daughter resonances. The helicity angle 
is unrestricted unless indicated otherwise. 



State 



inv. mass ( MeV) helicity angle 



K*j^+^-,K''j + ^^ 750 < ruKn < 1550 -0.85 < cose* < 1.0 



LO 



750 < < 1550 -0.80 < cos6' < 1.0 
470 < m^^ < 1070 -0.80 < cos0 < 0.80 
470 < m^^ < 1070 -0.70 < cos6i < 0.80 
735 < m-nTTTT < 825 

120 < rriTT, < 150 

488 < m^^ < 508 



asymmetry Ach = (r^ — r''')/(r^ + T^), where the su- 
perscript on the r corresponds to the charge of the B^ 
meson or the charge of the kaon for S° decays. 

The results presented here are obtained from data 
collected with the BaBAR detector [l^ at the PEP- 
II asymmetric-energy e+e~ collider located at SLAG. 
An integrated luminosity of 424 fb~^, corresponding to 
465 X 10^ BB pairs, was recorded at the T(4iS') resonance, 
with e+e^ ccnter-of-mass (GM) energy ^/s — 10.58 GeV. 

Gharged particles from the e^e~ interactions are de- 
tected, and their momenta measured, by five layers of 
double-sided silicon microstrip detectors surrounded by 
a 40-layer drift chamber, both operating in the 1.5-T 
magnetic field of a superconducting solenoid. We iden- 
tify photons and electrons using a GsI(Tl) electromag- 
netic calorimeter (EMG). Further charged particle iden- 
tification (PID) is provided by the average energy loss 
(dE / dx) in the tracking devices and by an internally re- 
flecting ring-imaging Gherenkov detector (DIRG) cover- 
ing the central region. 

We reconstruct i3-daughter candidates through their 
decays p° tt^tt^, /o(980) tt+tt", p+ 7r+7r°, 



K* 



^^°(iftt_n), K 



77, and 



K*° ^ K+n- 

TT+TT^. Gharge-conjugate decay modes are implied un- 
less specifically stated. Table U lists the requirements 
on the invariant masses of these final states. For the p, 
K* , and lo selections, these mass requirements include 
sidebands, as the mass values are treated as observables 
in the maximum-likelihood fit described below. For 
candidates we further require the three-dimensional flight 
distance from the primary vertex to be greater than three 
times its uncertainty. Daughters of p, K* , and a; candi- 
dates are rejected if their DIRG, dE/dx, and EMG PID 
signatures are highly consistent with protons or electrons; 
kaons must have a kaon signature while the pions must 
not. 

Table U also gives the restrictions on the K* and p 
helicity angle 6 imposed to avoid regions of large combi- 
natorial background from low- momentum particles. To 
calculate we take the angle relative to a specified axis: 
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for w, the normal to the decay plane; for p, the positively- 
charged daughter momentum; and for K* , the daughter 
kaon momentum. 

A i?-meson candidate is characterized kinemat- 
ically by the energy-substituted mass ttjes = 
\J (55 -I- Pq • Pb)'^/E^^ ~ pg and the energy differ- 
ence AE = Eg — ^\/s, where (£^o,Po) and {Eb,Pb) are 
four-momenta of the e+e~ CM and the B candidate, 
respectively, s is the square of the CM energy, and the 
asterisk denotes the e+e" CM frame. Signal events peak 
at zero for AE, and at the B mass for niES, with a 
resolution for AE (rriEs) of 30-45 MeV (3.0 MeV). We 
require \AE\ < 0.2 GeV and 5.25 < mes < 5.29 GeV. 

The angle 9t between the thrust axis of the B can- 
didate in the e+e~ CM frame and that of the charged 
tracks and neutral clusters in the rest of the event is 
used to reject the dominant continuum e^e~ ^ 99 (9 = 
u, d, s, c) background events. The distribution of | cos 9t\ 
is sharply peaked near 1.0 for combinations drawn from 
jet-like qq pairs, and is nearly uniform for the almost 
isotropic i3-meson decays. We reduce the sample sizes 
to 30000-65000 events by requiring Icos^tI < 0.7 for 
the ujp/fo modes and | cos^tI < 0.8 for the ujK* modes. 
Further discrimination from continuum is obtained with 
a Fisher discriminant ^ that combines four variables: the 
polar angles, with respect to the beam axis in the e~^e~' 
CM frame, of the B candidate momentum and of the B 
thrust axis; and the zeroth and second angular moments 
io,2 of the energy flow, excluding the B candidate, about 
the B thrust axis. The mean of is adjusted so that it 
is independent of the S-flavor tagging category 17| . The 
moments are defined by Lj ^J2iPi ^ jcos^ip , where 9i 
is the angle with respect to the B thrust axis of track or 
neutral cluster i and pi is its momentum. The average 
number of B candidates found per selected event in data 
is in the range 1.1 to 1.3, depending on the final state. 
We choose the candidate with the highest value of the 
probability for the B vertex fit. 

We obtain yields and values of Jl and Ach from ex- 
tended unbinned maximum-likelihood (ML) fits with in- 
put observables AE, toes> and, for the scalar, vector 
or tensor meson, the invariant mass and Ti. = cos 9. For 
each event i and hypothesis j (signal, qq background, BB 
background), we define the probability density function 
(PDF) with resulting likelihood £: 

V] = Vj{mEslV,{AE')Vj{r)Vj{m\,ml,Hm), (3) 
^ = ^^ET.Y^'P} ^ (4) 

i=l j 



m 



where Yj is the yield of events of hypothesis j, N is the 
number of events in the sample, and the subscript 1 (2) 
represents 37r (Kn or tttt). There are as many as three 
signal categories and the PDFs for each are split into two 
components: correctly reconstructed events and those 



where candidate particles are exchanged with a parti- 
cle from the rest of the event. The latter component is 
called self crossfeed (SXF) and its fractions are fixed to 
the values found in Monte Carlo (MC), (15-35)%. We 
find correlations among the observables to be small for 
qq background. 

From MC simulation 18] we form a sample of the most 
relevant charmless BB backgrounds (20-35 modes for 
each signal final state). We include a fixed yield (70- 
200 events, derived from MC with known or estimated 
branching fractions) for these in the fit described below. 
For i?+ Lup'^ we also introduce a component for non- 
resonant ujtt'^tt'^ background; for the other decays non- 
resonant backgrounds are smaller and are included in the 
charmless BB sample. The magnitude of the nonreso- 
nant component is fixed in each fit as determined from 
fits to regions of higher tttt or Ktt mass. For the cop 
modes, we also include a sample of 6 — s- c backgrounds; 
for the other modes, this component is not used since it 
is not clearly distinguishable from qq background. 

Signal is also simulated with MC; for the (i^7r)o line 
shape, we use a LASS model [3, [20| which consists of 
the i^^Q(1430) resonance together with an effective-range 
nonresonant component. For the /o(980), we use a Breit- 
Wigner shape with parameters taken from Ref. 2l|. 

The PDF for resonances in the signal takes the form 
Vi,sig{m\)V2,sig{m'^)QiH\,W2) with Q given by Eq. [1] 
or [21 modified to account for detector acceptance. For 
qq background we use for each resonance independently 
T^<i<i«.W = n9("^fc)n?(^fc): where Vg^iml) is a 
sum of true resonance and combinatorial mass terms. 
The PDFs for BB background have a similar form. 

For the signal, BB background, and nonresonant back- 
ground components we determine the PDF parameters 
from simulation. We study large data control sam- 
ples of D^Ti^ and _B+ D^p^ decays with 
I?° —> K^tt^tt'-' to check the simulated resolutions in AE 
and mESj and adjust the PDF parameters to account for 
small differences. For the continuum background we use 
(toeS; AE) sideband data to obtain initial values of the 
parameters, and leave them free to vary in the ML fit. 

The parameters that are allowed to vary in the fit in- 
clude the signal and qq background yields, Jl (for all 
VV and VT modes except 5° ^p'^), continuum back- 
ground PDF parameters, and, for ojp, the b ^ c back- 
ground yield. Since there is not a significant yield for 
^ bjp^ , we fix Jl to a value that is consistent with a 
priori expectations Q (see Table HI)) . For all modes ex- 
cept B^ —> ujp^ the signal and background charge asym- 
metries are free parameters in the fit. 

To describe the PDFs, we use simple functions such 
as the sum of two Gaussian distributions for many sig- 
nal components and the peaking parts of backgrounds, 
low-order polynomials to describe most background 
shapes, an asymmetric Gaussian for J-', and the function 
x'^ exp [—^(1 — a;^)] (with x 



mEs/Eg) for the 
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TABLE 11: Signal yield Y and its statistical uncertainty, bias Yo, detection efficiency e, daughter branching fraction product 
nSi, significance S (with systematic uncertainties included), measured branching fraction B with statistical and systematic 
errors, 90% C.L. upper limit (U.L.), measured or assumed /l, and Ach- In the case of wfo, the quoted branching fraction is a 
product with S(/o ^ tttt), which is not well known. {K-k)q refers to the S-wave Kn system. 



Mode 


Y 


Yo e 




B 


B U.L. 








(events) 


(events) (%) 


(%) (a) 


(io-«) 


(io-«) 







ujK*+ 

a; A"! (1430)+ 



101 ± 25 

72 ± 24 
8± 16 
540 ± 47 

191 ±36 
357 ± 39 
185 ± 32 

182 ± 30 
64 ±25 

30^ 
37^ 
411 ±43 



V21 
-18 
7+14 
-12 



O 1 A 

o ± 4 


1 r o 


oy.o 


A 1 


I.Z zh U.D ± U.Z 




n 70 _i_ n 1 /I -1- n no 
yj. ( Z dz U.i4dz U.Uz 


n A^. _Ln oc: i n no 
U.4o ± U.zo ± U.Uz 








2.0 


0/l_l_1 A_l_00 

Z.4 zh i.U =t U.z 




U.4i zh U.io =t U.UO 


n on -1- n qc: _l n no 
U.zy dz U.OO ztz U.Uz 


3 ± 2 


10.4 


29.7 


3.7 


4.8 ± 1.7 




0.37 ± 0.18 


0.22 ± 0.33 


0± 1 


13.6 


20.6 


0.5 


0.6 ± 1.2 




0.5 fixed 




49 ±25 


9.7 


59.5 


9.8 


18.4 ± 1.8 ± 1.7 






-0.07 ±0.09 ±0.02 








9.2 


27.5 ±3.0 ±2.6 






-0.10 ±0.09 ±0.02 


18 ±9 


6.4 


29.7 


5.9 


19.6 ±4.1 






-0.38 ±0.19 


34± 17 


9.1 


20.6 


10.6 


37.1 ±4.5 






-0.01 ±0.10 


19± 10 


11.9 


29.7 


5.0 


10.1 ±2.0 ±1.1 




0.45 ±0.12 ±0.02 


-0.37 ±0.17 ±0.02 








6.1 


21.5 ±3.6 ±2.4 




0.56 ±0.10 ±0.04 


0.14 ±0.15 ±0.02 


6±3 


8.2 


14.9 


7.2 


31.0 ±5.2 




0.52 ±0.10 


0.17 ±0.16 


10 ±5 


10.1 


10.3 


2.4 


11.2 ±4.9 




0.76 ± 0.26 


-0.04 ± 0.35 


-3 ±2 


9.5 


89.2 


1.9 


0.8 ±0.5 ±0.2 


1.6 


0.8 fixed 




1 ± 1 


14.4 


59.5 


4.5 


1.0 ±0.3 ±0.1 


1.5 






27± 14 


5.8 


89.2 


9.8 


15.9 ± 1.6 ± 1.4 




0.90 ±0.05 ±0.03 


-0.20 ±0.09 ±0.02 



Toes background distributions. These are illustrated for 
_B+ — > Lop'^ with projection plots of each fit variable in 
Figs. [U [2ji, and[3li. The parameters that determine the 
main features of the background PDF shapes are allowed 
to vary in the fit. 




g °0.74 0.76 0.78 0.80 0.82 8.0 0.2 0.4 0.6 0.8 1.0 



m,, (GeV) H, 



(e) 












-0.5 0.0 


0.5 







FIG. 1: Projections for B+ -> Lop+: (a) AS, (b) (c) 
msvr, (d) Ticj, and (e) 7i,r7r. Points with errors represent 
data and solid curves represent the full fit functions. Also 
shown are signal (blue dashed), b ^ c background (magenta 
dot-dashed), and total background (black long-dash-dotted). 
Charmless background and nonresonant background are too 
small to be seen. To suppress background, the plots are made 
with requirements on In £ that have an efficiency for signal of 
(40-60)% depending on the plot. 



We evaluate biases from our neglect of correlations 
among discriminating variables by fitting ensembles of 
simulated experiments. Each such experiment has the 
same number of events as the data for both background 
and signal; qq background events are generated from their 
PDFs while signal and BB background events are taken 
from fully simulated MC samples. Since events from the 
BB background samples are included in the ensembles, 
the bias includes the effect of these backgrounds. 

We compute the branching fraction B for each decay 
by subtracting the yield bias Yq from the measured yield, 
and dividing the result by the efficiency and the number 
of produced BB pairs. We assume that the branching 
fractions of the T(4S') to B^B^ and are each equal 

to 50%. In Table [ill we show for each decay mode the 
measured Jl, and Ach together with the quantities 
entering into these computations. For decays with K*'^ 
we combine the results from the two K* decay channels, 
by adding their values of — 21n£. For the significance S 
we use the difference between the value of —2 In £ for zero 
signal and the value at its minimum; the corresponding 
probability is interpreted with the number of degrees of 
freedom equal to two for modes with a measured and 
one for the others. For modes without a significant signal, 
we quote a 90% confidence level (C.L.) upper limit, taken 
to be the branching fraction below which lies 90% of the 
total of the likelihood integral in the region of positive 
branching fraction. In all of these calculations C{B) is 
a convolution of the function obtained from the fitter 
with a Gaussian function representing the correlated and 
uncorrelated systematic errors detailed below. 

We show in Fig. [5] the data and PDFs projected onto 
Toes- Figure [3] shows similar projections for the i^Tr and 
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niEs (GeV) 



25 5.26 5.27 5.28 5.29 

niEs (GeV) 
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(b)- 



















1.0 -0.5 0.0 0.5 



FIG. 4: B-candidate Ktv helicity projections for (a) loK*'^, 
and (b) uK*^ . The efficiency range and key for the curves 
are the same as for Fig. [2ja,b). 



FIG. 2; _B-candidate rriES projections for (a) ajK*", (b) loK""^ , 
(c) Lop^ /u>fo, (d) ojp^. The soHd curve is the fit function, black 
long-dash-dotted is the total background, and the blue dashed 
curve is the total signal contribution. For (a,b) we also show 
the signal components: K*{892) (red dashed), {Kt:)q (green 
dotted), and /('J (1430) (magenta dot-dashed). We show for 
(c,d) the fe — > c background (magenta dot-dashed), and for 
(c) the B° ujp° (red dashed) and B° cofo ( green dotted) 
components. The plots are made with a requirement on ln£ 
that has an efficiency of (40-60)% depending on the plot. 



TTTT masses. Figure S] gives projections onto Ti for the 
ujK* modes. 




08 1.0 1.2 1.4 

mK.(GeV) 



0.6 0.8 1.0 

m..(GeV) 



FIG. 3: B-candidate Ktt mass projections for (a) ujK*^, (c) 
LoK*'^, and tttt mass projections for (b) ujp^/ujfo, (d) ujp'^ . 
The efficiency range and description of the curves are the 
same as for Fig. [2] 



The systematic uncertainties on the branching frac- 
tions arising from lack of knowledge of the signal PDF 
parameters are estimated by varying these parameters 
within uncertainties obtained from the consistency of fits 



to MC and data control samples. The uncertainty in the 
yield bias correction is taken to be the quadratic sum 
of two terms: half the bias correction and the statistical 
uncertainty on the bias itself. We estimate the uncer- 
tainty from the modeling of the nonresonant and BB 
backgrounds by varying the background yields by their 
estimated uncertainties (from Ref. [16] and studies of our 
data). We vary the SXF fraction by its uncertainty; we 
find this to be 10% of its value, determined from stud- 
ies of the control samples. For the ifQ(1430) modes, we 
vary the LASS parameters within their measured uncer- 
tainties [l9[. For i?" ijjp^ where /l is fixed, the un- 
certainty due to the assumed value of /l is evaluated as 
the change in branching fraction when is varied by 
Ig'g. These additive systematic errors are dominant for 
all modes and are typically similar in size except for the 
error due to BB background, which is usually smaller 
than the others. 

Uncertainties in reconstruction efficiency, found from 
studies of data control samples, are 0.4%/track, 3.0%/7r'^, 
and 1.4%/ K';^ decay. We estimate the uncertainty in the 
number of B mesons to be 1.1%. Published data [lB| 
provide the uncertainties in the B-daughter branching 
fractions ( < 2%). The uncertainty in the efficiency of 
the cos^T requirement is (1.0-1.5)%. Since we do not 
account for interference among the K* components, we 
assign systematic uncertainties based on separate calcu- 
lations where we vary the phases between the three com- 
ponents over their full range. 

The systematic uncertainty on includes the ef- 
fects of fit bias, PDF-parameter variation, and BB and 
nonresonant backgrounds, all estimated with the same 
method as used for the yield uncertainties described 
above. From large inclusive kaon and B-decay samples, 
we estimate the Ach bias to be negligible for pions and 
—0.01 for kaons, due primarily to material interactions. 
Thus we correct the measured Ach for the K* modes by 
+0.01. The systematic uncertainty for Ach is estimated 
to be 0.02 due mainly to the uncertainty in this bias cor- 
rection. This estimate is supported by the fact that the 
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corrected background Ach is smaller than 0.015. 

In summary, we have searched for nine charmless 
hadronic S-meson decays as shown in Table [TTl and have 
observed most of them (for the first time in all cases 
except — > We calculate the branching frac- 

tions for wifg (1430) using the composition of {KTr)^ from 
Ref. We find B{B° luK^{U30)°) = (16.0 ± 

1.6 ± 1.5 ± 2.6) X 10-s and B{B+ ujKS{U30)+) = 
(24.0 ± 2.6 ± 2.2 ± 3.8) x 10"^ where the third er- 
rors arise from uncertainties in the branching fraction 
1^0(1430) ^ Ktt^ and the resonant fraction of (Kn)^. 
For most decays we measure Ach and find it to be consis- 
tent with zero. For VV and VT decays we also measure 
fr. For : Jl is near 1.0, as it is for B ^ pp 

jsj. For the VT B ^ ujK^{U30) decays /l is about 4cr 
from 1.0 for both charge states; it is similar to the value 
of ~0.5 found in B ^ <j>K* decays. Branching fraction 
results are in agreement with theoretical estimates Q ex- 
cept for the SV and VT decays where the estimates are 
more uncertain 11, 13l |. 
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